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Implicit Regularization (IReg) is a candidate to become an invariant framework in momentum 
space to perform Feynman diagram calculations to arbitrary loop order. In this work we present a 
systematic implementation of our method that automatically displays the terms to be subtracted 
by Bogoliubov's recursion formula. Therefore, we achieve a twofold objective: we show that the 
IReg program respects unitarity, locality and Lorentz invariance and we show that our method is 
consistent since we are able to display the divergent content of a multi-loop amplitude in a well 
defined set of basic divergent integrals in one loop momentum only which is the essence of IReg. 
Moreover, we conjecture that momentum routing invariance in the loops, which has been shown 
to be connected with gauge symmetry, is a fundamental symmetry of any Feynman diagram in a 
renormalizable quantum field theory. 
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I. INTRODUCTION 

A consistent renormalization program in QFT 
appeared after the work of Bogoliubov, Parasiuk, 
Hepp and Zimmerman (BPHZ) 0-0 in which 
a prescription to extract recursively the diver- 
gences of a multi-loop Feynman graph comply- 
ing with unitarity, locality and Lorentz invariance 
was presented. The BPHZ program generalized 
the Dyson's subtraction to general overlapping di- 
agrams to arbitrary loop order leading to the con- 
cept of renormalizable quantum field theory. Such 
program systematizes, according to the topology of 
the graph, the subtraction necessary to render the 
corresponding amplitude finite through the forest 
formula. The proof of finitude provided by the lat- 
ter is by construction regularization independent. 
However, for concrete predictions such as scatter- 
ing amplitudes in collision processes of elementary 
particles, the method of Dimensional Regulariza- 
tion (DReg) and minimal subtraction Q, [j| com- 
bined with Zimmerman's forest formula has proven 
to be an efficient and successful calculational tool 
particularly for gauge theories. The forest for- 
mula can be casted into a counterterm language 
by means of Bogoliubov's recursion formula [10| . 
complying with locality, Lorentz invariance, uni- 
tarity and causality. 
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To calculate S-matrix elements in a symmetry 
preserving fashion in a quantum field theoretical 
model sensitive to dimensional continuation on the 
space-time, the problem is more subtle. The con- 
struction of an invariant regularization framework 
is aesthetically more appealing but this is not the 
main motivation. Although in one hand impos- 
ing constraint equations derived from Ward iden- 
tities order by order in perturbation theory oblit- 
erates the need of an invariant regularization, on 
the other hand it renders the calculation more in- 
volved from the calculational viewpoint. Besides, 
if quantum symmetry breakings occur in pertur- 
bation theory, an invariant scheme is essential to 
judge it as physical or spurious. Supersymmetric 
gauge theories arc conspicuous examples of models 
in which regularization and renormalization play 
a fundamental role especially as new accurate ex- 
perimental evidence, viz. electroweak precision ob- 
servables [H[-[lI|, demands consistent theoretical 
calculations higher than one loop order to under- 
stand physics beyond the Standard Model. 

Therefore, the construction of an invariant reg- 
ularization is justified and, in order to be as reli- 
able as DReg (wherever DReg can be applied), it 
must be shown to comply with locality, Lorentz in- 
variance, unitarity and causality. Recently, an in- 
variant regularization framework (IReg) has been 
developed and shown to be consistent an d sy mme- 
try preserving in several instances [l5|-[3l|. The 
essence of the method is to write the divergences 
in terms of loop integrals in one internal momen- 
tum which do not need to be explicitly evaluated. 
Moreover it acts in the physical dimension of the 
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theory and gauge invariance is controlled by reg- 
ularization dependent surface terms which when 
set to zero define a constrained version of IRcg 
(CIReg) and deliver gauge invariant amplitudes 
automatically. Therefore it is in principle appli- 
cable to all physical relevant quantum field theo- 
ries, supcrsymmctric gauge theories included. A 
non trivial question is whether we can generalize 
this program to arbitrary loop order in consonance 
with locality, unitarity and Lorentz invariance, es- 
pecially when overlapping divergences occur. This 
is the main subject of this work in which we use 
the simplest renormalizablc field theoretical model 
to show how to implement IReg in such a way that 
it displays all the terms to be subtracted by Bo- 
goliubov's recursion formula automatically. All the 
other physical theories can be treated within the 
same strategy after space time and internal algebra 
are performed. Another result of this contribution 
is to show that if the surface terms are not set 
to zero they will contaminate the renormalization 
group coefficients. Thus, we are forced to adopt 
CIReg which is equivalent to demand momentum 
routing invariance in the loops. This feature leads 
us to conjecture that momentum routing invari- 
ance is a fundamental symmetry of any Feynman 
diagram. 



II. THE RULES OF IMPLICIT 
REGULARIZATION 

We restrict ourselves to massless theories and 
power counting infrared safe integrals. The first re- 
striction is justified because, as we show in Section 
IIV1 to implement a mass independent renormal- 
ization scheme in IReg we need only the massless 
basic divergent integrals that we present below. 
When infrared divergences do appear, a dual ver- 
sion of IRcg operating in coordinate space displays 
infrared divergences as basic divergent integrals as 
well, in a way that infrared and ultraviolet degrees 
of freedom are clearly distinguished [32^-[35j. 

Given the amplitude of a n-loop Feynman graph 
with L external legs, the basic strategy of IRcg is 
to free all divergences of external momenta and ex- 
press them in terms of basic divergent integrals in 
one loop momentum only. To achieve this purpose, 
we need to perform (n — 1) integrations, but the 
order in which they are performed is not clear a 
priori. In the next section we present a systematic 
way to choose the order of integration which, as 
a byproduct, displays the counterterms to be sub- 
tracted by Bogoliubov's recursion formula. Con- 



sidering that we made this choice, we can redefine 
the internal momenta in such a way that the inte- 
gral in ki is the 1-th we are going to deal with and 
it is typically of the form 

TVl ... Vm = f A^-^(h, qi ) / kf-^ \ 

J r»-%) 2 -M 2 ] I x 2 r 

ki 

(1) 

where I = 1 • • • n. In the above equation, qi is 
an element (or combination of elements) of the set 
{pi,...,p L ,k l+ i,...,k n }, J ki = fd d ki/(2ir) d and 
fj? is an infrared regulator. 

Since the original integral is infrared safe, the 
limit /i 2 — > is well-defined and must be taken 
in the end of the calculation. The logarithmical 
dependence appears because this is the character- 
istic behaviour of the finite part of massless ampli- 
tudes [HI . A is an arbitrary non- vanishing param- 
eter with dimension of mass which parametrizes 
the freedom one has to subtract the divergences 
(renormalization group scale). It appears at one 
loop level and survives to higher orders through a 
regularization independent mathematical identity 
(eq. [8|) as wc show in the end of this section. The 
function A Vl '" Vm (ki, qi) may contain constants and 
all possible combinations of fc; and qi compatible 
with the Lorentz structure. Care must be exer- 
cised when it contains a term like {ki—q.i) 2 . In this 
case, we must cancel it against one of the denom- 
inators because, as we are dealing with divergent 
integral s, sy mmetric integration is a forbidden op- 
eration [201, [37j . 

Now, we apply the rules of IReg. Assuming that 
a regulator A is implicit in the integral, we can use 
the following mathematical identity in the denom- 
inators: 

n (kl) -i 

1 = ' v (~iy(gf~2q l -k l y 

(h-q.f-y 2 ' ^ (fcf-^+l 

(ki) n (ki) 

(-l)" i ' ) ( g ?-2g i -fc t )"i' ) 
(A?-M 2 H*' ) [{h-qi) 2 -li 2 ] 

The values of are chosen such that all di- 
vergent integrals have a denominator free of qi. 

After the use of ([2]), the divergent integrals can 
be casted as a combination of 
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redefinition of a renormalization scheme. This fea- 
ture leads us to conjecture that momentum routing 
invariance is a fundamental symmetry of Feynman 
diagrams. 

At this point we notice that the divergences can 
be written in terms of one object namely 



In the above formula, the subscript log means 
that we are dealing with a logarithmic divergent 
integral (r = 2/3 — d). It is important to note that 
only this type of divergence appears because linear 
and quadratic divergent integrals vanish for mass- 
less theories [2{|, [381 ] . 

Although we have already reduced the diver- 
gences to basic divergent integrals free of external 
momenta, we can show that the integrals defined 
above are related. For example, in a case in which 
r = 2 we have 

CV> = £(§yf^{?4-V> 



2 J dkP 

k 



k" 



(fc 2 - y?yr- 



In 1 "' - 



k 2 - n 2 



(5) 



or equivalently, for short 



CV)-«r£«iC J ' +1 V) 



2V (Z-1)! 
d) 



(6) 



In the previous equation, T; is a surface term 
which is arbitrary and in general regularization de- 
pendent. It was shown in [lj|-[2(| that setting all 
surfaces terms to zero defines a constrained version 
of IRcg (CIReg) and corresponds to invoking mo- 
mentum routing invariance in the loops of a Feyn- 
man graph. This in turn is related to gauge invari- 
ance and it was shown that adopting CIReg is a 
sufficient condition to ensure gauge symmetry 31 1. 
One may verify that Dimensional Regularization 
(DRcg) evaluates the surface terms to zero what 
demonstrates that CIReg and DReg are compati- 
ble. In theories with less symmetry content such as 
scalar field theories one could ask whether momen- 
tum routing invariance plays any relevant role. We 
shall answer this question by calculating the first 
two coefficients of the </>g theory j3 function which 
are universal. We verify that the arbitrarity intro- 
duced by the surface terms cannot be hidden in the 
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(k 2 - n 2 ) d / 2 



hi 



i-i 



kf 



(7) 



However, the above integral is ultraviolet and 
infrared divergent as /i 2 — > 0. To separate these 
divergences and define a genuine ultraviolet diver- 
gent object we use the scale relation below [39[ 



'A\ ^ (-l) fc (/-!)! 



,kj^ ki (l-j) 

k=l v 7 j=l y J ' 



X 2 I ' 



(_l)d/2 



(47r)<V 2 r(d/2) 



(8) 
(9) 



Since we are dealing with infrared safe models 
the infrared divergence must disappear in the am- 
plitude as a whole. This in fact occurs because, as 
we use identity , the finite part of the amplitude 
will also have a logarithmical dependence in /x 2 and 
it is just the one expected to cancel the infrared di- 
vergence coming from the use of the scale relation. 
As mentioned before, we note that A parametrizes 
the freedom we have to subtract the divergences 
and becomes a natural candidate for a renormal- 
ization group scale. 

We repeat the above procedure until we are left 
with only one integral in the internal momentum 
and, consequently, we are able to express all diver- 
gences in terms of ^^ (A 2 ). One of the purposes 
of the next sections is to show, for a general n,-loop 
Feynman graph, how this program can be imple- 
mented in a systematic way which is compatible 
with Bogoliubov's recursion formula. 



III. SYSTEMATIC IMPLEMENTATION 
OF BOGOLIUBOV'S RECURSION 
FORMULA IN IREG 

In this section we develop an algorithm that im- 
plements IRcg to multi-loop Feynman graphs. It 
is constructed in such a way that it displays the 
terms to be subtracted by Bogoliubov's recursion 
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formula and thus fulfilling unitarity, locality and 
Lorentz invariance. 

In order to implement IReg in a systematic way 
to a n-loop Feynman graph we adapt identity @, 
which was initially conceived for one-loop order, to 
arbitrary order since it does not furnish us a nat- 
ural sequence in which the integrals must be per- 
formed. Therefore, our first task is to rewrite it in 
such a way that it evinces the divergent behaviour 
of the amplitude as the internal momenta go to in- 
finity in all possible ways. Restricting in ((5J) to 
be external momenta and expanding (pf — 2pi ■ k) 3 
with the familiar binomial formula yields, 

2(n (fc) -l) 

— L — £ / + (10) 



1=0 



where we defined, 

L//2J 



(k,Pi) - 



l-2j 



h { \ i ) W^WW-t ' 

(ii) 



/( , P „ 5 (-iH"(g-2 K . t) -r 

{k 2 ~p 2 )< [(fc — p^ 2 — ^i 2 ] 
6 ^ = \ 1 if x > ' 

B = n| fc ' + j — /, L :E J = niax{n G Z\n < x}. 

The terms /, p ^ are constructed in such a way 
that they behave like k~( l+2 ^ as k — >• oo and we 
choose the values of in order to assure the 
UV finitude of /( fe ' The above identity is the 
keystone of our procedure which, when applied to 
a given Feynman graph, can be summarized in the 
following steps: 

A. Identify the propagators which depend on the 
external momenta of the graph and apply iden- 
tity (HO]); 



B. Find out the minimum value of 



needed 



to assure the finitude of the terms that contain 
/ ( ki ' as ki — » oo in all possible ways; 

C. Repeat the above step for all propagators iden- 
tified in step [A] 

D. Identify the divergent terms and classify them 
according to all possible ways that the internal 
momenta approach infinity; 



E. Use the rules of IReg (in the way presented in 
Section [TTJ) in the terms identified in step [D] ac- 
cording to their classification; 

F. Set aside the divergent terms that contain 
/^(A 2 ) and apply the procedure again on the 
ones that do not. 

After stcpEl we have only two kind of terms: the 
ones in which ^g(A 2 ) multiplies an integral and 

the ones in which I® (A 2 ) multiplies only constants 
and/or polynomials in the external momenta. The 
first are just the terms cancelled by Bogoliubov's 
recursion formula while the latter are the typical 
divergence of the graph, i.e. after subtraction of 
subdivergences. In other words, our prescription 
implements IReg in a way that displays automati- 
cally the terms cancelled by Bogoliubov's recursion 
formula. 

We illustrate its applicability with some exam- 
ples. Since we are concerned only with the struc- 
ture of the divergences, we work with the simplest 
renormalizable quantum field theory: masslcss <f)\. 
In this theory, only graphs up to three external legs 
are divergent (Toj . The graphs with one external 
leg have only quadratic divergences and these al- 
ways vanish for massless theories (29|. Therefore, 
the graphs we deal with have only two or three ex- 
ternal legs and correspond to the renormalization 
of the propagator and the vertex functions respec- 
tively. 

In all the examples we present we choose a par- 
ticular momentum routing in order to simplify our 
calculation. However, we could have chosen a dif- 
ferent set of internal momenta and, if we have done 
so, we would obtain the same divergent structure 
expressed as basic divergent integrals. In other 
words, our prescription is not limited to a specific 
choice of momentum routing in the internal lines 
of a Feynman graph. 



A. One- and two-loop self-energy and vertex 
diagrams 

We begin with the one-loop correction for the 
propagator whose graph is 




FIG. 1. Graph P 
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The amplitude depicted by fig. Q] reads 



part of (fT5|) reads 



~ " 2 J (k 2 

k 



1 



1 



A* 2 ) [(& — J 3 ) 2 — M 2 



d 6 fc 
(2tt)8' 

(13) 



Notice that we have introduced an infrared reg- 
ulator in the denominators. We use identity f| 10[) 
in the denominator that contains the external mo- 
mentum to obtain 



3(1) 

~~F 



(fc 2 -/i 2 ) 



2(n (k) -l) 

E • 

1=0 



(k,p) 



■ (fc, p) 



(14) 



We must choose rv- k > in order to guarantee the 
finitude of the term that contains By 
power counting we find that rS k ^ > 2 and thus 

(fc, p) 



we choose rv- k ' = 3. Having found the value of n 



we can extract the divergent terms. As /; goes 
like fc _ (' +2 ) we find by power counting that they 
are given by: 



1. Quadratic divergence 



f 1 
Jo 



(fc, P) 



(fc 2 - )1 2 ) J (fc 2 - fl 2 ) 



(15) 



2. Linear divergence 

(fc, p) 



J(k 2 -^ 2 ) J] 



2p ■ fc 



(fc 2 -fi 2 ) 3 



(16) 



3. Logarithmic divergence 

■(k, P ) 



(fc 2 - fi 2 ) J (fc 2 - n 2 ) 3 

k k 



(2p-k) 2 
(fc 2 - II 2 ) 



(17) 



The quadratic and linear divergences vanish in 
the limit /i 2 — » 0. The logarithmic one contains 
two terms which can be identified with Ij^i^ 2 ) 
and Ii og (pL 2 ) respectively. We use identity Jfj|) to 
express all divergences in terms of // og (/i 2 ) and 
therefore the amplitude is given by 



"(1) „2 

^ = -^i log (» 2 ) + 2 P 2 r 1 

9 6 



finite, 



(18) 



where Ti is an arbitrary constant steaming from a 
surface term. The explicit expression for the finite 



(fc 2 - ^ 2 ) 



2 7 (fc 2 - /x 2 ) 4 
fc 

= ^lnf- r 
6 



9 



(fc — p) 2 — /i 2 
0(/i 2 ). (19) 



Using the scale relation (eq. [8]) and taking the 
limit ji 2 — > wc finally obtain 



2 2 

7 (i)__sV 



7io 9 (A")-66lii(-^]+-^-12Ti 



(20) 




In a similar way we obtain the amplitude of the 
one-loop correction for the vertex function: 



Iio 9 (X 2 ) - &e In(-fi) + 26 6 - h(pi,pa) 

FIG. 2. Graph V (1) 

In the amplitude depicted by V^, h(pi,p2) is a 
function of pi and p 2 which vanishes if P2 = 0. 

Although the previous examples are too simple 
to display the terms to be subtracted by Bogoli- 
ubov's recursion formula (the one-loop graphs do 
not contain subdivcrgences) , we presented them 
here to familiarize the reader with the rules of IRcg 
and identity (fTO)) . We proceed now to the two- loop 
corrections for the propagator and the vertex func- 
tions. The graphs needed in the renormalization of 
the propagator are: 





FIG. 3. Graphs PJ ; and P£' respectively 

?( 2 ) w 



The amplitude corresponding to P\ is given by 

tt( 2 ) 1 r 

^- = - A(fc 1 )A(fc 1 -p)A(fc 1 -fc 2 )x 
W 2 J 



k\kn 



A(fc 2 )A(fc 2 -p), (21) 



G 



where we defined 

A(fc 4 



(22) 



We begin by applying identity (fTH)l in the prop- 
agators that depend on the external momenta to 
obtain 



/A(fci)A(fci -fc 2 )A(fc 2 )x 



~ 2 ( n ( fc i>-l) 



2(n< fc 2)_i) 



m=0 



(23) 



Our next task is to determine the values of n^^. 
They are chosen in order to assure the fmitude of 
the terms that contain /( fc i>p) as fc, — > oo in all 
possible ways. We consider first The terms 

that contain p) are 



k±k 2 



A(fc 1 )A(fc 1 - fc 2 )A(fc 2 )/~ 



"2(n( fe 2'-l) 
m— 



(fcl, p) , 



"(*2,p) 



A(fc 1 )A(fc 1 -fc 2 )A(fc 2 )/^>rfA(fc 2 -p). (24) 

fclfe2 

We want to guarantee that the above integral is 
finite as fci — > oo. There are two cases: 

1. Finitudc as fci — s- oo and fc 2 fixed: vS kl ' > 0, 

2. Finitudc as k\ — > oo and fc 2 — > oo: n' fcl ' >2, 

which leads us to conclude that rS kl ^ should be at 
least 3. In a similar fashion, we obtain n^ 02 * 1 = 3. 

Having found the values of rS ki \ we proceed to 
identify the divergent terms contained in (|23jl as 
fci and/or fc 2 go to infinity. There are three pos- 
sibilities. We start with the case fci — >• oo and fc 2 
fixed where the divergence terms are of the type 



A(fci)A(fc 2 )A(fci - fe)/ 



(fci,p). 



(k 2 ,p) 



+ f 



0=2, p) 



(25) 



As f ' kl ' pl ) goes like k^ l+2 \ we find by power 
counting that the divergent terms in this case are 
given by 



Af = / A(fci)A(fc 2 )A(fci -fc 2 )/ 



(fei.p), 



j (*2, p) _|_ j (k 2l p) 

_m=0 

= y A 2 (fci)A(fci -fc 2 )A(fc 2 )A(fc 2 -p). (26) 

/cifc 2 

We consider now the case where fc 2 — > oo and fci 
is fixed. Repeating the previous reasoning, we find 
that the divergent terms are 



Af= I A(fci)A(fc 2 )A(fci-fc 2 )/ (fc2 ' p) : 



^ /; ( fe ,P) + /(fcl ,p) = 

j=o 

/" A 2 (fc 2 )A(fci-fc 2 )A(fci)A(fci-p). (27) 



A:i k 2 



Finally we consider fci — > oo and fc 2 — > oo si- 
multaneously. The choice of nA^'J = 3 (i = 1,2) 
assures us that the divergent terms must be of the 
type 



A(fci)A(fc 2 )A(fci - k 2 )f l {kl ' p) f^ p) . (28) 



k± k 2 

By power counting, we obtain that I and m arc 
constrained by / + m < 2. The cases I = and 
in = 0,1,2 are contained in Af (eq. [26]) while the 
cases m = and Z = 0,1,2 are contained in A^ (eq. 
[27]) . We are therefore left with the case I = m = 1 
which reads 

Af =/ A(fc 2 )A(fci -fc 2 )A(fci)/ 1 (fcl ' p) /i (fc2 ' p) = 

fcl &2 

= y A 3 (fci)A(fci -fc 2 )A 3 (fc 2 )(2p-fci)(2p-fc 2 ). 

(29) 

Summarizing, the divergent terms are: 
1. Divergences as fci — > oo and fc 2 is fixed 

Af= yA 2 (fci)A(fci-fc 2 )A(fc 2 )A(fc 2 -p), (30) 

k± Ai2 
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2. Divergences as fc 2 — > oo and fci is fixed 
Af= J A 2 (k 2 )A(k 1 ~k 2 )A(k 1 )A(k 1 -p) 1 (31) 

fclfe2 

3. Divergences as fci — > oo and fc 2 — > oo simul- 
taneously 

Af = [ A 3 (fc 1 )A(fc 1 -fc 2 )A 3 (fc 2 )(2p-fci)(2p-fc 2 ). 



(32) 



Therefore, the divergent content of i s given 
by Af + A 2 + Af — Af . The last term corresponds 



to the case (I 



0) 



(fci, p) n (k 2 , P) 



Af= I A(fc 2 )A(fci-fc 2 )A(fci)/ ^/ ' 
= J A 2 (k 2 )A(k 1 -k 2 )A 2 (k 1 ) 

ki k2 



(33) 



and must be subtracted because it is counted twice. 

The above classification of the divergent terms 
in different cases (the term Af can be thought of 
as the intersection between the cases k\ — >• 00 and 
fc 2 fixed, fc 2 — > 00 and fci fixed) furnishes us a 
natural order in which the integrals must be per- 
formed and allow us to implement IReg to multi- 
loop Feynman graphs in a systematic way. More 
importantly, as a byproduct it also displays the 
terms to be subtracted by Bogoliubov's recursion 
formula as we shall verify. Roughly speaking, for 
each of the cases we have studied characterizing 
the divergent behaviour in all possible ways that 
the internal momenta go to infinity we may read- 
ily apply the formalism developed in Section HTl to 
define basic divergent integrals. For this purpose 
we use identity <j2j) to each case taking fc; to be the 
internal momentum that we pick to go to infinity. 
In our example, fci( 2 ) in and k\ and fc 2 in 

Examining Af and Af , we notice that both have 
the same structure and the integral in which we are 
going to use the rules of IReg is given by 



A 2 (fc,)A(fc, - kj), i,j = 1,2 and i 56 j (34) 



This is the same amplitude of graph V\ (fig H]) if 
we identify p\ — > kj and set p 2 = 0. Thus we can 



readily write 
Af = Af + 



af, i,j = 1,2 and i^j 
Af = J A{k 3 )A{k 3 -p) [l log {\ 2 )}, 



A(fc,)A(fc J - p) In (-%J^-\ . (35) 



We turn to Af ■ Since the integral in fci is finite, 
we evaluate it by Feynman parametrization. We 
insert the result in the integral in fc 2 and use the 
rules of IReg to obtain 



a 3 = A t = b &P 2 
Similarly 



WA 2 ) 



- j A^(fc 2 ) 

fc 2 




2Ti 



H-n 2 
A 2 



(36) 
(37) 



in the limit p 2 —> 0. 

We will show that the terms Af (i = 1,2) are 
just the ones which must be subtracted by Bogoli- 
ubov's recursion formula. Let us set them aside 
for the time being and evaluate the rest (af). Us- 
ing identity (|T0|) in the propagator that depend on 
the external momentum, we identify the divergent 
terms. In the limit /x 2 — > 0, the only one that 
contributes is 



A(fc,)/, 



(kj , p) 



-be In - 



26 6 
(38) 



which can be expressed by 



af = b 6 p 2 



& 2 ) 



-J ios (A 2 ) + 8Ti-4T 2 



(39) 



12) 

Hence, the divergent content of phis surface 
terms is given by 



?(2)c 



19 4 



(af + a 2 + af + Af + A% t 



(40) 



The two last terms are just the ones to be sub- 
tracted by Bogoliubov's recursion formula. In fact, 
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in order to subtract the subdivergences of this par- 
ticular graph we must add the following countert- 
erms 





• = (— 1) x Divergence of /s. 



(->) 



FIG. 4. Counterterms for P 



whose amplitudes are, respectively 
^ [ A(fc 2 )A(fc 2 -p) [-WA 2 )] = -£ (-Af] 



(41) 



^ J A(fcx)A(* x -p) [-^ os (A 2 )] = (-If) . 
fei 

(42) 

Notice that we are adopting a "MS" scheme that, 
in IRcg, corresponds to the subtraction of basic 
divergent integrals [i~8l |. 

Therefore, subtracting the subdivergences yields 



n( 2 ) 

_ _ 

«5 4 " 6 



2lS(A a )-^J lofl (A a ) 
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54Ti - 24T 2 + finite 



(43) 



We turn now to the two loop nested graph (Pg ) 
whose amplitude is 



"(2) 



1 



z.g 4 2 



A 2 (fci)A(fc x -p)A(k 2 )A(k 1 -k 2 ). 



(44) 



and fci fixed does. They are given by 



A 2 (fc!)A(fc! -k 2 )A(k 2 ) 



fci fc 2 



P)_j_ y (fci, p) 



A 2 (/ci)A(fc! - fc 2 )A(fe 2 )A(fci -p). 



(45) 



fcifcg 



For definiteness call the above integral Bf . This 
is the only one which we have to deal with (the di- 
vergent terms from the case k\ — > oo and k 2 — > oo 
simultaneously are contained in the above inte- 
gral). One may notice that it is just the original 
amplitude of the graph but now we have a natural 
order to implement IRcg. We use its rules in the 
integral in k 2 to obtain 



Bf = Bf+f3f 



Bf = I A(fci)A(fci-p) 



3 v ; 



ff = J A(ft 1 )A(* 1 -p) 

k 



be, ( kj-jl 



+ J A(*i)A(Ai -p) 



4Ti 



(46) 



We apply the procedure again in j3f to find the 
following divergent terms 



Pf=J A(fcx) 

fci 



(fe, P) 



.1=0 



A(fci) 



LZ=0 



(fe 2 ,p) 



106 f 



4Ti 



7 Zo9 (A 2 )-46 6 T 2 



4Ti 



2& fi 



(47) 



We use identity (|10| in the propagator that de- 
pend on the external momentum and find that the 
choice n( fcl ) = 3 guarantees the finitude of the 
terms that contain /( fel >p) as k\ — > oo in all pos- 
sible ways. We proceed to identify the divergent 
terms. The case k\ — » oo and k 2 fixed does not 
contain any divergent term while the case k 2 — > oo 



Notice that an arbitrary valued surface term 
appears multiplied by a divergence expressed by 
hog{^ 2 )- Should we have adopted CIReg which 
sets T; = as required by momentum routing in- 
variance (and gauge symmetry) it would not have 
appeared. We shall keep all surface terms until 
the end to see whether they play any role in the 



9 



physics of a less symmetrical theory such as scalar 
field theories. 

Therefore, the divergent content of zr B ' plus sur- 
face terms is given by 



1 



— (2)oo 



(48) 



Again, the last term is just the one which is 
subtracted by Bogoliubov's recursion formula since 
the counterterm we must add is 



i£/A(fa)A(Ai-p)[iJ loi ,(A a )] 

fci 



fci 



(— l)x Divergence of 



fci 



fci 



FIG. 5. Counterterm for P, 



(2) 



«.9 4 



After we add the counterterm we obtain 
b 6 p 2 



6 



^o 3 (A 2 ) + ^-3T! 



+ finite | . 
(49) 



Notice that the term proportional to 4Ti/j og (A 2 ) 
is not subtracted as a subdivergence since we are 
adopting a "MS" scheme in IReg. In the next sec- 
tion we discuss what would happen if we had in- 
cluded the surface term in the counterterm. 

At this point, we can write down the renormal- 
ization of the propagator at two loop order 

"(2) _ ~(2) , ~(2) 
- — "A T 



AP J 56 6,(2) ri 2N 296 6 



*3 S — ^ (A') ~ (A 2 ) - 20T 2 - 



-2Ti 



2^og (A 2 



776 fi 



6Ti 



finite j . 
(50) 



We now turn to the renormalization of the ver- 
tex. The graphs we need to evaluate are: 





fci fc'2 

Pi P2 pi k 2 p2 pi P2 

FIG. 6. Graphs \ Vj, 2) and V^ 2) respectively 



(2) 

The amplitude depicted by V\ is given by 



A 



(2) 



e y > A(fc 1 )A(fc 2 - fci) JjA(*j - pi)A(^ - pa). 



fci fc 2 



(51) 

We apply identity (fTQ|) in the propagators that 
depend on the external momenta and make the 
choice n^ 1 ^ = r4 = n^ 2 ' = n^ 2 ' = 1 in the 
usual way to obtain 

A(fc 1 )A(fc 2 -fc 1 )[/ (fc - Pl) +/( fc -^) _ 



fci fc'2 



f( kl ' P2 ~> -j- f ( fc i> Pa) 

J^(fc2, P2) _|_ J (fc 2 , p 2 ) 



jj^ 2 ' pi) _|_ J (fc 2 , pi) 



(52) 



The divergent terms comes only from the case 
fc 2 — > oo and fci fixed and they are given by 

A$= J A(h)A(A 2 -fc 1 )/ (fa ' Pl V (fa,, * ) >< 

fci fc'2 



(fci. Pi) _|_ j (fci, pi) 



jjfci, Pa) _|_ j (fci, p 2 ) 



=y > A(fc 1 )A(fc 2 -fc 1 )A 2 (fc 2 )A(fc 1 - Pl )A(fc 1 -p 2 ) 

A: i A.: 2 

= A£ + at 

A-i = y > A(fc 1 )A(fc 1 -pi)A(fcj - pa) [/, og (A 2 )] , 
fci 

e 2 A 2 



a£ =^ A(fcx) A(fcx -pOA(*i -p 2 ) In (- 
fci 



(53) 

We apply the procedure again in to yield the 
divergent term below 

o^ee fe 6 /A(fc 1 )/ (fc - Pl) / (fcl ' P2) f2-ln('- fc? ~^ 



fci 



A 2 



26 6 // os (A 2 )-6 6 /g(A 2 ) 



(54) 
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Hence, the divergent content of is given by 



Repeating the procedure in yields 



AW°° = -ig 5 [at + At] (55) 
where the last term is cancelled by the counterterm 




-ig 5 J A(fci)A(fci -pi)A(fc a -p 2 ) [- W A ')] 



Pi Vl 



= ( — l)x Divergence of 




pi — /ci pi — k\ 



(2) 



FIG. 7. Counterterm for 

Therefore, after the subtraction of the subdiver- 
gcncc we have 



A 



(2) 



-ig 5 



( A 2 ) + 2I log (\')+ finite . (56) 



(2) 

We turn to graph V B whose amplitude is 



A 



(2) 

B — 



-ig* 



i J A 2 (k 1 )A(k 1 - Pl )A(k 1 -p 2 )x 

x A(fc 2 )A(fc 2 - fci). (57) 



fcifc 2 



We use identity (fTU|) and choose n^ 1 — n 2 = 1- 
The divergent terms are all contained in the case 
k 2 — > oo and fcj fixed and they are given by 

Bf = y > A 2 (fc 1 )A(fc 2 )A(fc 2 - fc x )x 
fcife 2 

=jA 2 {k 1 )A{k 2 )A(k 2 -k 1 )A{k 1 - Pl )A{k 1 - P2 ) 



5 A = J A(fci)A(h -pi)A(*i -p 2 ) 
fei 

=y A(fc!)A(fci - Pl )A(*i -p 2 )x 



— (A ) 



fci 



4T 1 + |ln 



A 2 



86e 
9 



(58) 



pi = f A(*i)n/ ' 

I i=l 



(fcl,pi) 



In 



A 2 



4Ti 

9 



= y^S(A 2 )+4T 1 / io9 (A 2 ) -^ os (A 2 ) (59) 

and, therefore, the divergent content of A^ plus 
surface terms is given by 



A 



(2)oo 



+B$). 



(60) 



The last term is subtracted by Bogoliubov's 
recursion formula since the counterterm for this 
graph is 



/ A(fci)A(fci -pi)A(fci -pa) [| W A ')] 




• = (— l)x Divergence of 



FIG. 8. Counterterm for 



Subtracting the subdivergence we have 



A 



(2) 



—E=J I io a i X >-{-^ 



2T 1 )I log (X 2 )+ finite. 

(61) 



(2) 

Finally, we are going to evaluate the graph VA . 
Calling its amplitude A c we obtain 



A 



(2) 

C — 



= i jA(k 1 )A(k 2 )A(k 1 -p 1 )A(k 2 -p 2 )x 

k%k 2 

A(fci + k 2 - px)A{ki + k 2 - pa). 

(62) 



As usual, we choose = n 2 k ^ = n[ kl+k2 * > = 
^(fci+fc 2 ) _ ^ an( j g n( j on jy divergent term 

comes from the case k± — > oo and fc 2 - >■ oo simul- 
taneously and it is given by 

/~tA _ / r (ki, pi) f (fe 2 , pi) r (fci+fc 2 , pi) ,. (ki+k 2 , p 2 ) 

b i — / Jo Jo Jo Jo 



ki k 2 



A 2 (fc 1 )A 2 (fc 2 )A 2 (fc 1 + fc 2 ). 



(63) 



fclfc2 
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After the integration over k 2 and the use of the 
rules of IReg we have 



A ( V = A<? = -tg% [l log {\ 2 



finitel 



(64) 



Collecting all the results we find that the renor- 
malization of the vertex at two loop order is given 

by 

A' 2 ' eeA ( , 2) +A (2) -a(2) 



A 



5fo 6 T (2) n 2 
'Y 1 log\ A 



+ ig 5 x (finite) 



6T1 Ilog{\ 



(65) 



Although our procedure found success in all the 
examples already presented, one may wonder how 
general it is. To answer this question, wc apply it 
to diagrams with more than two-loops. 



B. General Algorithm 

In order to see the generality of our method we 
consider the graph below 




whose amplitude is given by 



^(4) 

-ig 8 



A^ftiJA^fa - P )A{k 2 )A{k 2 -k 1 )x 

k i 

A(k 2 -k 3 )A(k 3 )A(k 3 -h)x 
A(fc 4 )A(fc 4 - fci +p). (66) 



As usual, we use identity (jT0|) in the propaga- 
tors which depend on the external momenta and 
choose n( fcl ) = rS ki ~ kl ^ = 3. Next, we must iden- 
tify the divergent terms as the internal momenta 
go to infinity in all possible ways. For this particu- 
lar graph, we are lead to an ambiguity because the 
cases fc 2 (3) ~~ ^ 00 an d ^1, &3(2)) ^4 fixed; k 2 , k 3 — > oo 
and k\ , fc 4 fixed contain the same divergent terms 
which read 

r 4 i 2 

A 2 (fci) 



' 4 
.1=0 



j (fci, V) _|_ J (fci, p) 



A(fc 2 )x 
k 3 )A(k 3 )A(k 3 - ki)x 



A(fc 2 -fci)A(fc 2 

" 4 

A(fc 4 ) ^/ r (^- fe i>p) 



r=0 



(67) 



Therefore, it is not clear if we must sum all of 
them (and subtract the terms counted twice) or if 
we must consider just one of the cases. However, it 
is not the main problem. Supposing we agree what 
case(s) we must consider, our next task would be 
to apply the rules of IReg (as presented in Section 
UlJ) according to the previous classification. In any 
of the cases, the integrals in fc 2 ( 3 ) are performed 
and we are unable to identify the counterterm 



= -if- J A 2 (fci)A(fci-p)A(fc 2 )A(fc 3 )x 
A(fc 2 -fci)A(fc 2 - fe 3 )A(fe 3 - fci) 




Hence, we need to develop a more general proce- 
dure which will contain steps lA"llFl as a subcase. We 
present this new procedure (which can be casted as 
an algorithm) and, afterwards, explain it step-by- 
step through the previous example. The algorithm 
reads 



YES- 



NO 

6 



r 

{5} YES • 



a 



YES NO 



YES 
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NO 



YES 



YES 



YES ■ 



NO 



y + i 



NO 



NO 
_L_ 



x+1, y 
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■NO 
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1. Let S e be the set of subgraphs Gi that share 
one external leg with the whole graph G. Let 
S e be the complementary set. Is S e an empty 
set? 

2. Apply the procedure stated before (steps 
EHF} to the amplitude. 

3. Identify the propagators which depend on 
the external momenta and do not belong to a 
subgraph contained in S e . Use identity (fTOjl 
in such propagators and find out the values 
of rij in the usual way. 

4. Set x = and y = 1. Group all internal mo- 
menta in the set £^=0 an< ^ define its com- 
plement (which is an empty set) by K^=$ ■ 



the ones in which Ii og {\ 2 ) multiplies an integral 
(that correspond to the terms cancelled by Bogoli- 
ubov's recursion formula) and the ones in which 

Iio g {\ 2 ) multiplies only constants and/or polyno- 
mials in the external momenta (that correspond to 
the typical divergence of the graph). 

As stated, we are going to clarify the method 
with the previous example. In the following, we 
present the result we obtain after each one of the 
steps. 

• Step : answer is NO 




5. Is there a subgraph Gi whose internal mo- 
menta are precisely the ones grouped in 



6. Are there divergent terms as all internal mo- 

- (y) 

mcnta contained in JCx go to infinity and all 
elements of are kept fixed? 

7. Have all the divergent terms identified in step 
[5] been stored? 

8. Is K-x a subset of a already stored K^f (i = 

0,1,.-. ,3- 1)? 

9. Store the terms identified in step [5] as well as 

(if fC^ is a subset of a already stored 
, erase the result corresponding to the 
latter) . 

10. Choose a new combination of x numbers of 
internal momenta and group them in the set 
K, x . Is it possible? 

11. Choose x numbers of internal momenta and 
group them in the set /ci y ' ) . Is it possible? 

When the algorithm gets to asterisk (*) , we 
write down all the results stored. Each one is 
classified according to a particular set of internal 
momenta that go to infinity. If this set has two 
or more elements, we identify the subgraph that 
contains them and run the algorithm on it. On 
the other hand, if the set has only one element 
we use the rules of IReg (as presented in Section 
HTj) . At this point, we set aside the terms that 
contain ifjgi}?) and apply the algorithm on the 
others. Eventually, we obtain two kind of terms: 




} 




s e = 



• Step H 
n( fcl >=3, A(fci-p) 

• Step [Ml 



U=0 



K 



(i) 



{&!■■■ fc 4 }, K\ 



(1) 



• Step []||] : answer is YES 

The internal momenta k\ • ■ • fej (which are con- 
tained in IC^Sq^) belong to a subgraph of the set 



• Step M : answer is YES 
The divergent terms are 



A 2 (fci) 



k\ •••ht 



'<2 



(fell P) f Ol, p) 



A(fc 2 )x 



E ftfi 
i, v= o 

A(fc 2 -fc 1 )A(fc 2 - fc 3 )A(fc 3 )A(fc 3 - h)x 
A(fc 4 )A(fc 4 -fa+p) 
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Step [jfj] : answer is NO 
Step M 



We store the divergent terms given by 



well as K%-q ] 



• Step |[T0"1| : answer is NO 

• Step [jnJ] : answer is YES 
We choose ki and obtain 

• Step [0] : answer is NO 

None subgraph contains {k 2 ■ • ■ fc 4 } as its inter- 
nal momenta. 



Step [TO] : answer is YES 



We choose k 2 and obtain 

ic^ ] = {h , * 3 , fc 4 }, /ciri 1} - {fc 2 } 

• Step [0] : answer is NO 

None subgraph contains {fci, fc 3 , fc 4 } as its inter- 
nal momenta.. 



Step [10] : answer is YES 



We choose fc 3 and obtain 

tt 1 ' = {h , *2, fc 4 }, /cir^ - {fc 3 } 

• Step [5] : answer is NO 



None subgraph contains {fc 4 , k 2 , fc 4 } as its inter- 
nal momenta. 



Step [TO] : answer is YES 



We choose fc 4 and obtain 

K^ 1 1) ={h---k 3 }, = {**} 

• Step [0] : answer is NO 

None subgraph contains {fci ■ ■ ■ fc 3 } as its inter- 
nal momenta. 



Step QJJ] : answer is NO 



Step [jTT]l : answer is YES 



We choose k\ , k 2 and obtain 



K. 



(y=i) 



{fc 3 ,fc 4 }, K 



(y=i) 



{ki,k 2 } 



• Step [5] : answer is NO 



None subgraph contains {fc 3 ,fc 4 } as its internal 
momenta. 



Step [TO] : answer is YES 



We choose ki , fc 3 and obtain 

ft 1 ' = {k 2 ,h}, 4r 2 x) = {fci, a 3 } 

• Step [M] : answer is NO 

None subgraph contains {fc 2 ,fc 4 } as its internal 
momenta. 



Step [TO] : answer is YES 



We choose ki , fc 4 and obtain 

= {fc 2 , fc 3 }, /cir 2 x) = {Ai, a 4 } 

• Step []H]] : answer is YES 

The second subgraph of S e has {A2,A3} as its 
internal momenta. 

• Step [||] : answer is YES 
The divergent terms are 



A 2 (fci) 



f( k i, p) + j (fci, p) 



Z=0 



A(fc 2 )x 



A(k 2 -h)A(k 2 - fc 3 )A(fc 3 )A(fc 3 - Ai)x 
A(fc 4 )A(fc 4 - fci (69) 

• Step [0] : answer is NO 

One may notice that (|69|) possesses more terms 
than (1681). 



StepH 
ince 

result corresponding to the latter. Therefore, we 



Since /Q,^_ 2 is a subset of K^Lq , we erase the 



14 



store only the divergent terms given by as well The divergent terms are 



as K.^ =2 ' = { k 2,h}- 



0=1) 

x=2 

• Step [[TU1| : answer is YES 
We choose k 2 , fc 3 and obtain 
(v= 2 ) -su. u,\ rO=i) 



• Step O : answer is NO 

None subgraph contains {k\,k 4 } as its internal 
momenta. 

• Step [Ml : answer is YES 



We choose k 2 , &4 and obtain 

£f= 2) = {ki, k 3 }, = {k 2 , k 4 } 

• Step H_J : answer is NO 

None subgraph contains {fci,fc 3 } as its internal 
momenta. 



Step [10] : answer is YES 



We choose fc 3 , k 4 and obtain 



0=2) 



{fci,fc 2 }, /c£l- 2 2) = {fc 3 ,fc 4 } 



Step [M] : answer is NO 



None subgraph contains {k\,k 2 } as its internal 
momenta. 



Step [10] : answer is NO 



Step [UJ : answer is YES 



We choose k\,k 2 , fc 3 and obtain 



K, 



0=i) 

x— 3 



= {fc 4 }, /Cits 15 = {fc 1 ,fc a ,* 9 } 



• Step [0] : answer is YES 

The first subgraph of S e has {fc 4 } as its internal 
momentum. 

• Step : answer is YES 



y p) _(_ j ( k ^p) 



2=0 



A(fc 2 )x 



A(k 2 -h)A(k 2 - fc 3 )A(fc 3 )A(fc 3 - fci)x 

A(fc4)A(fc4-fc!+p). (70) 

• Step [M] : answer is YES 



One may notice that (f70|) contains the same 
terms of P?]). 

• Step ffl : answer is NO 

Since we erased the result corresponding to 
K^Tq , we notice that K^=3 is not a subset of 
a previously stored K^=P ■ 

• Step m 

We store the divergent terms given by ([TOf as 



well as /C 



0=1) 

a:— 3 



Step [M| : answer is YES 



We choose ki,k 2 , k 4 and obtain 

4r 3 2) = {fc 3 }, /cir 3 2) = {fci,fc2,M 

• Step [||]] : answer is YES 

The fourth subgraph of S e has {fc 3 } as its inter- 
nal momentum. 

• Step [||] : answer is YES 
The divergent terms are 



-1 2 



.1=0 



A(fc 2 )x 



A(fc 2 -fc!)A(fc 2 - fc 3 )A(fc 3 )A(fc 3 - fcOx 

A(fc4)A(fc4-fc!+p). (71) 

• Step H : answer is YES 

One may notice that ([7T|) contains the same 
terms of tiffih and d70l). 



Step [Ml : answer is YES 
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One may notice that k4= 3 is a subset of ■ 



Step [JO] : answer is YES 



We choose fci, fc 3 , fc 4 and obtain 

& = {fc 2 }, 4r 3 2) = {fci,fc 3 ,fc4} 



Step [M] : answer is YES 



The third subgraph of S e has {fc 2 } as its internal 
momentum. 

• Step : answer is YES 

The divergent terms are 



~i 2 



i=0 



A(fc 2 ): 



A(fc 2 -fc 1 )A(fc 2 - fc 3 )A(fc 3 )A(fc 3 - fci)x 
A(fc 4 )A(fc 4 -fc 1 +p). (72) 

• Step [M] : answer is YES 

One may notice that (|72|) contains the same 
terms of (IS9| and (l70l) . 



• Step [M] : answer is YES 

One may notice that K^Z^ is a subset of K^S^ ■ 



Step [10] : answer is YES 



We choose fc 2 , fc 3 , &4 and obtain 



fe 2) ={fcl}, 4l- 3 2) ={fc 2 ,fc3,fc4} 



• Step [M] : answer is YES 

None subgraph has {fci} as its internal momen- 
tum. 



Step [10] : answer is NO 



Step [TTJ : answer is YES 



We choose k\, fc 2 , fc 3 , k 4 and obtain 

- 0, /dfjf - {Ai, fe, * 3l M 



• Step [5] : answer is NO 



Step [10] : answer is NO 



Step [TTJ : answer is NO 



At this point, we write down the results stored: 



!• K-x=2 ] = {>2,fc 3 } and cq. ([69]) 



A[ 4) = J A 2 (k 1 )A 2 {k 1 -p)A(fc 4 )A(fc 4 - ki + p)x 



x A(fc 2 )A(fc 2 -fc 1 )A(A :2 ~fc 3 )A(A :3 )A(A :3 -A :i ), 



2. /C 



(y=i) 



- x =3 - {M and eq. ([70]) 
A 2 4) = y A 2 (fci)A 2 (fc! -p)A(fc 4 )A(fc 4 -ki+p)x 



x A(k 2 )A(k2-k 1 )A(k 2 -k 3 )A(k 3 )A(k 3 -k 1 ), 

3. Term counted twice (intersection between 
the cases UJ and ©) 

4 4) = / A 2 (fci)A 2 (fc! -p)A(k A )A(k 4 -k 1 +p)x 



ki ■■■kA 



x A(k 2 )A(k2-k 1 )A{k 2 -k 3 )A(k 3 )A(k 3 -k 1 ), 

The term is divergent as two momenta 

(k 2 ,k 3 ) go to infinity. Therefore, we identify the 
subgraph that contains these momenta and apply 
the algorithm on it. The result can be written us- 
ing a graphic notation as below (ioj 



A(fc 2 )A(fc 2 -fc 1 )A(fc 2 -A :3 )A(fc 3 )A(fc 3 -fc 1 



= ig 



1 kfbe 
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2^(A')-f // os (A 2 ) + finite 



3 

J2 I AikJAiki-h) 



hog{\ 



2\ 1 



Therefore, can be re-expressed as 



16 




+ 




We set aside the terms that contain /^(A 2 ) and 
run the algorithm again on the others. The results 







We turn to the term A^\ Since it is divergent 
as only one momentum goes to infinity (Au), we 
just apply the rules of IReg to obtain 

^ jA(k 4 )A(k i -p') = -^p^-[l log (\ 2 ) + finite] 

&4 



where p' = k\ — p. 

Thus, A% is given by 




We are left with the term A^ which we can be 



re-expressed as 




or equivalently 




+ 





+ 

P V 






= £ a„( P )/^(A 2 ) + finite 



+ — F- 



A brief comment: although we have not cal- 
culated the finite part of P A (which appears in 

A^), it can be shown to be of the form logarithm 
plus constants, allowing us to express all diver- 
gences in terms of (A 2 ). 

At this point, we notice that all divergent terms 



contain 1^ (A 2 ) and the divergent content of 5 
is given by 

„(4)oo _ ig" \ A (4) 



^(4) 



' 4 



4 4) - 4 4) 






As stated before, the first term (in which -%1(A 2 ) 



multiplies only constants and/or polynomials in 
the external momentum) corresponds to the typ- 
ical divergence of the graph while the others (in 
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which //^(A 2 ) multiplies an integral) correspond 
to the terms to be subtracted by Bogoliubov's re- 
cursion formula. 



C. iV-loop vertex and self-energy diagrams of 
the ladder type 



To conclude this section, we focus on n-loop 
graphs. We applied our method in two examples: 
the first contained only nested divergent subdia- 
grams while the other contained all kinds of sub- 
divergences including overlapping ones. In both 
cases, our method found success since it was able to 
display the terms to be subtracted by Bogoliubov's 
recursion formula and the typical divergence of the 
graph was also written in terms of a well defined set 
of basic divergent integrals in one loop momentum 
only. Due to the lack of space, we just present the 
result of the first graph while the other one is dis- 
cussed in more detail. The nested diagram reads 




FIG. 9. V (n) 



and its amplitude, after the subtraction of the sub- 
divergences, is given by 



A(") 



+ l„2n+l — m log 



r^ff (A 2 ) + nnite, (73) 



tained 

(0) 



1, 



°m=l — x ! 
An-1) 



L m^l — 



2 ^E 



p—m—l 



\ Fp-m+l ( n -2) 



D n = 



F„ 



p — m + 1 



1 



e=0 



(2-e)(l-e) 
1 

(e + 2)(e+l)(e-l) 



(74) 



Before we proceed to the other graph, we state 
an important result which was essential to show 
that the typical divergence of is expressed in 
terms of BDI's (the procedure to compute finite 
integrals can be found in |4lj ) 



jA(k 3 -k a )A 2 (k 3 ) 



E c ir 1)ln * 



pi—i 
i+i 



E c ir 1)J ff( A2 ) + E c £ v 



A 2 



/ k 2 -// 5 
1 K a A* 



m— 1 



(75) 



Now we are ready to evaluate the graph with 
overlapped subdivergences which reads 




FIG. 10. Graph P (n) 
Its amplitude is given by 



- n n 

A / J]A(fc r )A(fc r -p)n A ( fc ;-^-i), 

(76) 



1 = 2 



where A 



where the coefficients Cm can be recursively ob- graphic notation suitable for our purposess 



Before applying our method, we introduce a 
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ki--k 



E cjr 1) in m - 1 



m — 1 
n— j+1 



E 



( 71 -J')l„m-1 

77. m 



A 2 (k i )A 2 (k j )A(k j - 

n A(/c i )A(fc ; -fc, 

I=i+1 



-4 J 



E 



m — 1 
n— J+1 



(™-j)] n m-l 



E > 

777 — 1 

A 2 (fc J )n 1 A(fc i )A(fc i+1 



-4 J 

ki — k 



E 



, (i-l) l n rrL-l 



777—1 

n— j+1 



E cir-^in— 1 



777,-1 



A 2 (fc) II A(fc ; )A(fc, 



E 



n— j+1 

E - 



, l n m-l 



j n m-l 



fcf — M' 



k 2 -fi* 



kj-i) x 
i)A(*,- 



P) 



fc^ — u' 
_J _ 

— 



fc ; )A(fc 



i)A(fc 



kt—p' 



k 2 - —ii' 

3 _ 





-P) 



P) 



= A 



n A(fc r )A(fc r -p) n A(fc i -fc i _ 1 

r — i l=i-\-l 



(-1) x E c^I&'CA 3 ) 



777,-1 



(-1) x E ^-^/^(A 2 ' 



n / A(fc r )A(fc r - P ) n A(fc ! -fc ! _ 



(-i)x e^Oa 2 ; 



n— j + 1 

E • 



(n-i)] n m-l f k j 



IT 



fl A(fc r )A(fc r -p) n A{h-h-x) 

r=i+l ;=i+2 




= -4 / 



(-1) x E c£- 1J OA a 

m — 1 

E c^ln™" 1 



n A(fc r )A(fc r - P ) n A(fe-fc,_i) 

r—i £=i+l 



and define 



ki—k 



if 



A 2 



A 3 (fc,)x 



A 3 (fc,)x 



n-j+l / u2_.2 

E ^In-M-^ 

m— 1 \ 



l=i+l 



(77) 



In all the terms defined above, we consider that 
a bad-defined product is valued to unity. We are 
now in a position to use our algorithm but, since all 
divergent subgraphs share one leg with one of the 
external legs of the entire graph, it reduces to steps 
DOfI Therefore, we just apply identity (fT0|) in the 
propagators that depend on the external momenta, 
choose nS ki ' = 3 (i = 1 • ■ • n) and find the following 
divergent terms: 

1. Divergences as k\ — > oo and k<i ■ ■ ■ k n arc 
fixed 





+ 



2. Divergences as k n — > oo and ki ■ ■ • k n —i arc 
fixed 



- 1 




nY- = - - l 




+ - 1 



71-1- 



3. Term counted twice (intersection between 
the cases HJ and (J2J) 
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<|>-<|p-<[> 



+ 




+ 



n-1 



To prove that the above statement is correct, we 
use mathematical induction. We begin by applying 
our method to a typical term of the first summa- 
tion and obtain the following divergent terms: 



1. Divergences as k t 
are fixed 



oo and k 



i — a+1 



4. Divergences as ki ■ ■ ■ k n — y oo simultaneously 

I 1, n 

At this point, we would set aside the terms that 
contain Ii og (X 2 ) and use the method again in the 
others. Although this is the basic procedure to 
be followed, adopting it will prove to be a harder 
way to find the typical divergence of the graph for 
arbitrary n. Therefore, instead of applying our 
method in a blind way we may notice that the 
results above were obtained after the use of (|75|l 
which means that we performed some of the n in- 
tegrals. In the first two cases we performed just 
one integral (hi and k n respectively), while in the 
third one we performed the integrals in k\ and k n . 
Hence, if we had performed only one integration, 
the divergent content of would be given by 




+ 




+ 



n-1 



Now, we suppose that performing i — 1 integrals 
lead us to 



i-1 /T\ i-1 
J2 - \i-a\ \n-a\ - + £ (-1) 
o=0 vJJT a=l 




i k-1 

X/ £ Ik—a,n—a+l 
k=2 0=1 



+ 6(i-2)£ (-l)x 
r=2 o=l 




+ 



i-1 

£(-i)x 

a=l 




+ 










i — a+1 




n — a 







2. Divergences as k n - a — y oo and ki 
are fixed 



kji — a- 





a+1 







- i — a 




n — a— 1 







3. Term counted twice (intersection between 
the cases HJ and ©) 




4. Divergences as ki- a ■ ■ ■ k„- a — > oo simultane- 
ously 

1 i—a, n — a 

In other words, the divergent content of a typical 
term of the first summation is given by 




n — a - + - i — a+1 



n — a - — - 1 — a 




'a + 1 



+ -i — a 




n — a — I- - -(i — a) (n — a) h „_„ 



FIG. 12. 

We consider now a typical term contained in the 
second summation. As we seek the divergent con- 
tent of P( n ) after i integrals are performed, we are 
able to rc-cxpress it as below 





i — a+1 





+ 



i — a+1 



FIG. 11. 



FIG. 13. 
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Inserting the results of (fig [T2"|) and (fig [T5|) at 
(figllip we obtain 



E 

£J=0 



+ E (-i)x 

a=l 




i+1 k-1 

+ EE/ A;— a, n— a+1 + 

fc=2 a=l 



©(< - 1) E E (-i) x 

r=2 o=l 




+ 





E(-i) : 

a=l 



One may notice that the result above is just 
the one we supposed (fig QT]) if we set i — > i + 1 . 
Therefore, by mathematical induction, we prove 
our statement and we will be able to identify the 
typical divergence of the graph. Performing n — 1 
integrals lead us to 



n-l 

E 

a=0 




n-l 

+ E (-i) x 

a=l 



n k—1 

+ EE/ k— a, n— a+1 + 

k=2 a=l 



n—1 r—1 
+ 0(n-2)£ £ (-1) x 
r=2 a=l 




+ 



n-l 

E (-i) x 

a=l 




+ 




give a null contribution. The third summation can 
be evaluated and the result is 



n k — 1 k—1 a 



k—1 a—1 m—1 m' — l 
m-\-m — 1/ 1 \ / 



""^ 71 V (m + m'-2)! ( TO+m '-Q a2 >, 
(ro + m'-l-i)! J ^ lA j 



g — "-fc ) 

d « = 2b 6 d 

k—q 



d" 



d(ea)5 



(i-i) A - 1 



1 



Gk-q, 



(79) 



The others terms are just the ones that are sub- 
tracted by Bogoliubov's recursion formula since 
the counterterms for the graph we are dealing with 
can be re-expressed as 







We apply our method in the first summation and 
find that the divergent terms coming from it are 
given by 



n—l n— a a+1 

f 2 EE E 

a— m—1 m' — l 
M-l 



, (n-a-1) (a) 



(M-l) ^ 2 
lay 



(A 2 ) 



M = m + m'. 



(78) 



In the second summation all the terms contain 
only quadratic divergences and, therefore, they 




a loops 



(— l)x Divergence of 



Therefore, the typical divergence of PW (after 
the subtraction of subdivcrgenccs) is given by the 
sum of with (|79]) . 

One may notice that our method once again dis- 
played the terms to be subtracted by Bogoliubov's 
recursion formula and allowed us to express the 
typical divergence of the graph in terms of well 
defined BDFs. Therefore, we claim that the pro- 
cedure we developed here (which is summarized in 
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the algorithm) is the systematization of IReg to 
multi-loop Feynman graphs. Although the results 
we already obtained are restricted to massless the- 
ories, we can develop a similar procedure to treat 
massive ones as shown in the next section. 



IV. MASSIVE THEORIES 

We generalize identity (|10|) to massive theories 
and implement a mass independent regularization 
scheme. The procedure we develop here allows us 
to write the divergences of any theory (massive or 
not) in terms of the set of basic divergent integrals 
we have already presented and, therefore, justifies 
our restriction to massless theories in the previous 
sections. 

The identity we are going to use is 



1 



(k—pi) 2 — m 2 — [i 2 
where we defined. 



hl k ' Pi) +h^' Pi \ (80) 



1=0 



[1/2] 
3=0 



1-3 



{-p 2 +m 2 y{2 Pi -k) 

(k 2 - p2)J+l-j 



pQ = (-!)"• \pj-m 2 -2 Pi -kY 



(81) 
(82) 

(fc 2 - A i 2 )«r'[(fc-^)2_ m 2_ (U 2] 

Notice that, if we set m 2 — > 0, we recover iden- 
tity (fTO)) . The procedure to be followed is similar 



(fc, Pi) 
l 



was 



to the one of the massless theory since h 
also constructed in such a way that it goes like 
k~( l+2 > as k — > oo. The main difference is that we 
must apply identity ([50]) in all propagators which 
depend on the external momenta and/or the mass. 
In order to familiarize the reader with this identity, 
we evaluate the one-loop contribution for the prop- 
agator 



| / A m (fc)A m (fc-p) 




FIG. 14. Graph M (1) 

In the above picture we defined (42 

1 



A m (fc) 



k 2 — m 2 — /i 2 



(83) 



We apply (|50"|) in all denominators to obtain (we 
defined po = and p\ = p) 



2(n< fc) -l) 

E 

1=0 



Ak, 0) 



h ^ °) 



2(»i 



(fe) 



E 

q=0 



(84) 



As usual, we choose the minimum value of n\ 
(i = 0, 1) in order to guarantee the finitude of the 
terms that contain h ( k ' as k — > oo. In this case, 



we find that n 



(k) 



(k) 



3. At this point we 



identify the divergent terms: 
1. Quadratic divergence 



h (k, o), (fe, 



2. Linear divergence 



h (k,o) h (k, P ) 



k k 
3. Logarithmic divergence 



{k 2 -/i 2 ) 2 ' 



2p-k 
{k 2 - n 2 f ' 



h?' 0) h o h ' p) 



(k 2 - n 2 ) 3 ' 



(85) 



(86) 



(87) 



(2p-k) 2 (p 2 -m 2 ) 



(k 2 -n 2 ) 4 {k 2 -n 2 ) 3 



One may notice that the quadratic and linear 
divergences are the same of the massless case and 
they give null contributions. Therefore, after the 
use of the rules of IReg, the remaining terms yield 



<r 



WA 2 ) + yTi 



finite . 
(89) 

The two-loop corrections are obtained in the 
same way and we just present the results. Aftcr 
the subtraction of the subdivergences, the two-loop 
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overlapped graph {P\ ' ) furnishes 



ig 4 



©A 2 ) 



-^// os (A 2 ) + 9T 1 -4T 2 



hp 2 



2ll 2 >(X 2 ) + 5I loa (X 2 ) 



log 



bam 2 



finite, 



(90) 
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(2) 

while for the two-loop nested graph (Pg ') we ob- 
b 6 p 2 b 6 ', 



tain 




~(2m) 






{(- 


ig* 






fbb 6 


+ 


^27 




2b 6 T 2 


•( 


3 



18 



2TA 2 



l4 2 i(A 2 



2& 6 _ T \ 2 

— +6Tijm 



^ os (A 2 ) 



+ 2(Ti) 2 p 2 + finite 



(91) 



Therefore, the renormalization of the propagator 
at two-loop order is given by 



~(2m) 



5b 6P 2 



56 fi m 2 



\ 18 2 

^296 6 2TA o 



4 2 ](a 2 ) 



I 54 3 / V 3 
106 6 T 2 776 6 T! 



+6Tibi' 



WA 2 ) 



9 



2(Tj 



finite. 



(92) 



To conclude, we calculate the renormalization 
group functions. As usual, we make the definitions 



Zmm, g a = Z g g, 



(93) 



which allows us to divide the bare lagrangian in 
two parts: the renormalized lagrangian and the 
counterterms. Explicitly, 



L = \ [(<^) 2 - 



V- 



+ 1 -[A(d^r-Bm^}+^, 
A = Za — 1, B = Zd,Z m — 1, C = Z„Z 



(94) 



We also define the renormalization group func- 
tions by 

A<91nZ dg x dlnZ g 

^=2~dX' /3 " A 9A = -- 9A ^A-' (95) 

7™ = --A^=A^. (96) 
m aA aA 

Supposing that A, £? and C have a expansion in 
the coupling constant g as below 

A = A ig 2 + A 2 g 4 + 0(g 6 ), (97) 
B = B ig 2 + B 2 g A + 0(g 6 ), (98) 
C = C l9 2 + C 2 g 4 + 0(g% (99) 



we obtain 



In Z<f, =Aig 2 



A 2 - 



)ff 4 + 0( ff 6 ), (100) 



lnZ m =(B 1 -A 1 )g 



A 2 + #-^V 4 + 0(.9 6 ), (101) 



lnZ„ = 



S 2 - 



d - — ) g 



Since 7, 7 m and /3 can also be written as 



7 = 7i.9 2 +72<? 4 + 0(ff 6 ), 



7m = 7^V+7W + 0(5 6 ), 

/3 = /3ig 3 + /V + o( 5 6 ), 

we deduce the relations below 

- X dAl 

71 ~ 2 9A ' 

A3 r a 2 

72 = ^ 1+ 2aA( v A2 -f 
7 W =A A( J B 1 -A 1 ), 



(103) 
(104) 
(105) 



y (2) 



dD 

2(B 1 -A 1 )(3 1 +\—, 

Ds ( B2 _ A2 + |_|), 



i? = Co - 



3A 2 3A 2 C 2 



(106) 
(107) 

(108) 
(109) 
(110) 
(111) 
(112) 
(113) 



Because we chose a "MS" scheme in IReg, the 
coefficients Aj, £>i and Cj are given by (equations 
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(EH}, (E2D, §B and fig. HD 

Ai = -\ho g {\ 2 ), 

Bl = -lhog{\ 2 ), 



(114) 

9 (A 2 ), (H5) 
(116) 



B, = -^S(A 2 ) + + eT. ^^A 2 ), (117) 



Ci = ^-MA 2 ), 
5& 6 r (2) 

2 



(118) 

C = - 5 4C(* 2 ) + (^ + eA^A 2 ). (119) 



of the f3 function which are known to be univer- 
sal. Hence, we must set them to zero even in a 
scalar theory. This in turn corresponds to assure 
momentum routing invariance in the loops of an 
arbitrary Feynman diagram [171 ] and allow us to to 
conjecture that momentum routing invariance is a 
fundamental symmetry of Feynman graphs. 

It is also important to note that when anoma- 
lies come into play, they may express themselves, 
from the perturbative viewpoint, as an explicit de- 
pendence in momentum routing as pointed out by 
Jackiw [43j . In such case, to exhibit democratically 
the anomaly among the Ward identities, surface 
terms should be left arbitrary 
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Since we have the relations 

di log (\ 2 ) _ 2x2 di log (\ 2 ) 



X- 



dX 



dX 2 



-266, (120) 



A- 



<|(A 2 ) 



OX 

we finally obtain 



= -2I log (X 2 ) - 36 6 , (121) 



7 = 



13ff 4 ig 4 



12(4tt) 3 432(4tt) 6 3(4tt) 3 



(122) 

5g 2 97 g 4 16»g 4 6 

7m " 6(4^)3 + 108(4tt)6 + 3(4tt)3 1 + U(9 h 

(123) 

3o 3 125g 5 5ig 5 „ fis 



4(4tt) 3 144(4tt) 6 (4tt) 3 



(124) 



One may notice the appearance of arbitrary sur- 
face terms in the renormalization group functions. 
At first glance, this feature may be a result of the 
definition of the "MS" scheme in IReg which cor- 
responds to the subtraction of basic divergent in- 
tegrals only. However, if we use a different scheme 
in which we subtract basic divergent integrals and 
surface terms we obtain the same result above for 
the first two coefficients of the (3 function. 

As pointed out earlier, surface terms are related 
to gauge and supersymmetry in such a manner 
that setting them to zero guarantees the invari- 
ance of the amplitude regarding these symmetries. 
In our present case we are not dealing with gauge 
or supersymmetric theories but with a scalar one. 
Therefore, it is natural to ask if the surface terms 
will play any role and the answer is that they in- 
troduce an arbitrariness in the first two coefficients 



V. CONCLUDING REMARKS 

We have shown that IReg being a strong candi- 
date for a symmetry preserving invariant regular- 
ization is consistent to locality, Lorentz invariance, 
unitarity and causality. This was achieved through 
an algorithm that implements IReg to multi-loop 
Feynman graphs in such a way that the terms to 
be subtracted by Bogoliubov's recursion formula 
are displayed automatically. We have also demon- 
strated that CIRcg (which corresponds to setting 
the surface terms to zero and is the sufficient condi- 
tion to deliver gauge and supersymmetric invariant 
Green' s functions) should be adopted in theories 
with less symmetry content as well. We learn from 
this that momentum routing invariance is a funda- 
mental symmetry of any Feynman diagram and all 
regularization procedures should comply with it at 
the expense of bringing arbitrary non physical pa- 
rameters into the amplitude. In our example this 
is manifest in the two-loop coefficient of the /? func- 
tion. It is important to notice that IReg is a n-loop 
invariant program that displays the divergences as 
basic divergent integrals. If one evaluate the BDFs 
or surface terms using a specific regularization (di- 
mensional, Pauli-Villars, etc) then it becomes ap- 
parent that any regularization which attributes a 
non zero value to surface terms could crash with 
momentum routing invariance and thus assign a 
non physical value to the universal coefficients of 
the P function. An exceptional case corresponds 
to anomalies in perturbation theory which mani- 
fests themselves as a breaking of momentum rout- 
ing invariance [43| . Then, for instance, in order to 
democratically display the anomaly between the 
vector and axial sectors of the AW triangle (AB J 
anomaly) surface terms must be let arbitrary as 
physical free parameters. 
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